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OO , evidence that similar approaches are relevant to A/ = 8 supergravity, we investigate 

^ ' possible iterative structures and applications of Wilson loop techniques in maximal 

supergravity. We first consider the two-loop, four-point MHV scattering amplitude 
^ ■ in A/" = 8 supergravity, confirming that the infrared divergent parts exponentiate, 

and we give the explicit expression which represents the failure for this to occur for 
the finite part. We observe that each term in the expansion of the one- and two- loop 
amplitudes in the dimensional regularisation parameter e has a uniform degree of 
transcendentality. We then turn to consider Wilson loops in supergravity, showing 
that a natural definition of the loop, involving the Christoffel connection, fails to 
reproduce the one-loop amplitude. An alternative expression, which involves the 
metric explicitly, is shown to have a close relationship with the physical amplitude. 
We find that in a gauge in which the cusp diagrams vanish, the remaining diagrams for 
this Wilson loop correctly generate the full one-loop, four-point J\f = 8 supergravity 
amplitude. 



^ {a. brandhuber, p. j. heslop, a. nasti, w.j. spence, g.travaglinijOqmul. ac.uk 



1 Introduction 



Evidence of recursive structures in the S-matrix of gauge theories has emerged in the 
past few years. In 2003 Anastasiou, Bern, Dixon and Kosower (ABDK) [1] made 
the remarkable observation that the planar, four-point MHV scattering amplitude in 
jV=4 supersymmetric Yang-Mills theory at two loops can be written as a polynomial 
of the one-loop amplitude, plus a kinematic-independent numerical constant. Subse- 
quently, innovations prompted by twistor string theory, and in particular improved 
generalised unitarity techniques, have made possible the calculation of higher-loop 
amplitudes both in Yang-Mills and in gravity. Bern, Dixon and Smirnov (BDS) [2] 
were able to show that the iterative structure uncovered in [1] holds up to three loops, 
and put forward a conjecture for the all-loop, ra-point MHV amplitude in A/" = 4 su- 
per Yang-Mills at the planar level, in which the all-loop amplitude is obtained by 
a suitable exponential of the one-loop amplitude multiplied by the cusp anomalous 
dimension. 

Iterative structures were first discovered by analysing the soft and collinear be- 
haviour of amplitudes in gauge theory [3-10]; the remarkable fact uncovered in [1,2] 
was that the finite parts of the MHV amplitude also follow the same pattern induced 
by the expected exponentiation of the infrared divergences. The BDS proposal was 
checked at three loops in the four-point case in [2], and subsequently in [11] for the 
two-loop, five-point amplitude. In a very recent paper [12], a discrepancy was found 
between the form of the amplitude conjectured by BDS and an explicit two-loop cal- 
culation of the six-point amplitude. The result at six points shows that the structure 
is that of a polynomial in the one-loop amplitude, plus a kinematic-dependent finite 
remainder function. 

In a related development [13], Alday and Maldacena proved the correctness of the 
BDS proposal for the four-point amplitude at strong coupling using the AdS/CFT 
correspondence0 In their calculation, the exponentiation of the one-loop amplitude 
occurs through a saddle point approximation of the string path integral, which in 
the AdS case turns out to be exact. Furthermore, they showed that the computation 
of amphtudes at strong coupling is dual to the problem of finding the area of a 
string ending on a lightlike polygonal loop embedded in the boundary of AdS space. 
This, in turn, is equivalent to the method for computing a lightlike polygonal Wilson 
loop at strong coupling using the AdS/CFT correspondence, where the edges of the 
polygon are determined by the momenta of the scattered particles. In a subsequent 
paper [15] the same authors showed that the BDS conjecture should be violated for 
a sufficiently large number of scattered particles. Further evidence of a breakdown of 
the BDS conjecture was also found in [16]. 

-'^For a recent review, see [14]. 
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The work of [13] suggested that the calculation of a Wilson loop with the same 
polygonal contour could be related to that of the MHV scattering amplitude even at 
weak coupling. This was proved in [17] for the one-loop four-point A/" = 4 amplitude, 
and by three of us in [18] for the infinite sequence of one-loop MHV amplitudes in 
J\f = 4 super Yang- Mills. This surprising Wilson loops/amplitudes duality was later 
confirmed at two loops for the four- [19], five- [20], and six-point case [21,22]. On 
the Wilson loop side, exponentiation naturally emerges as a result of the maximal 
non-Abelian exponentiation theorem [23,24]. Furthermore, the form of the four- 
and five-point expression of the Wilson loop is determined (up to a constant) by an 
anomalous dual conformal Ward identity [20] , and found to be of the form predicted 
by the BDS ansatz. A similar dual conformal symmetry was found for the integral 
functions appearing in the expression of the multi-loop amplitudes in [25]. Since 
conformal invariance is not restrictive enough to fully constrain the n-side polygonal 
Wilson loop for n > 6, it was perhaps not surprising that at precisely six points the 
BDS conjecture turned out to be incorrect [12]. It is intriguing however that the 
Wilson loops/amplitudes duality does not seem to break down - indeed, the results 
of [12] and [22] show numerical agreement between the Wilson loop and the six-point 
gluon amplitude at two loops. 

These iterative structures in gauge theory and string theory have been found at 
the planar level. Planarity appears to be a key ingredient of the story - for instance, 
the non-planar parts of the four-point MHV amplitude at two loops do not respect 
the same iterative structure as the planar part [1]. Planarity would also appear to 
be an important ingredient in any relation to integrability - the cusp anomalous 
dimension appearing in the BDS proposal is also determined by an integral equation 
derived in [27] using integrability. An analytical solution to this equation was recently 
presented in [28]. 

It is natural to ask if gravity shares any of these remarkable properties. Gravity 
is a non-planar theory, hence it is perhaps even more unexpected to find regularities 
in the higher-loop structure of its S-matrix. However, the mounting evidence of 
interconnections between the maximally supersymmetric theories of A/" = 4 Yang- 
Mills and J\f = 8 supergravit}|§ gives reason to be more optimistic. Perhaps the 
potentially most impressive similarity between these two theories is the conjecture 
that the J\f = 8 theory could be ultraviolet finite [30-36], just like its non-gravitational 
maximally supersymmetric cousin. Furthermore, gravity is also well understood in 
the infrared thanks to the results of [37], where it was found that infrared singularities 
can be resummed to the exponential of the one-loop infrared divergences, in complete 
similarity to those of QED [38,39]. 

With these motivations in mind, in this paper we would like to initiate a twofold 
^The paper [29] reviews the subject up to 2002. 
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investigation in A/" = 8 supergravity. Our first goal will consist in looking for possible 
iterative structures and cross-order relations using the known results at one and two 
loops for the MHV four-point scattering amplitudes. We confirm that the infrared- 
divergent parts exponentiate, but we observe a failure for this to occur for the finite 
parts, in contradistinction with the four- and five-point amplitudes in A/" = 4 Yang- 
Mills. On the other hand, we find that, similarly to the A/" = 4 MHV amplitude, 
each term in the expansion of the one- and two-loop A/" = 8 MHV amplitudes in the 
dimensional regularisation parameter e has a uniform degree of transcendentality (or 
polylogarithmic weight). This is very intriguing, and leads to the speculation that 
maximal transcendentality [40] could be yet another common feature of A/" = 4 super 
Yang- Mills and J\f — S supergravity. 

Our second aim is to investigate possible relationships between gravitational scat- 
tering amplitudes and gravitational Wilson loops. This second objective is further 
motivated by some calculations of gravity amplitudes in the eikonal approxima- 
tion [41,42], and by our belief that there should exist a strong link between the 
eikonal approximation [43-45] (performed in specific kinematic regions) and the more 
recent polygonal Wilson loop calculations (performed without reference to any specific 
kinematic region). 

The rest of the paper is organised as follows. In the next section we will describe 
the known one- and two-loop MHV amplitudes in A/" = 8 supergravity, and use them 
to show that the two-loop amplitude, minus one half of the square of the one loop 
amplitude, is finite, consistently with general arguments concerning the exponentia- 
tion of infrared divergences in gravity. We give the explicit expression for this finite 
term. 

In Section 3 we turn to a one-loop Wilson loop calculation. One candidate for the 
Wilson loop expression, given by an integral of an exponential involving the Christof- 
fel connection, is shown not to give the one-loop supergravity amplitude correctly. A 
second expression for the gravity Wilson loop is then studied, motivated by its ap- 
plication in the eikonal approximation to gravity. This involves the metric explicitly 
and is not gauge invariant, however the failure of gauge invariance is restricted to 
terms localised at the cusps of the Wilson loop. 

The individual cusp diagrams and finite diagrams have the structure expected for 
the A/" = 8 MHV amplitude (with the tree-level amplitude stripped off); however, 
after summing over all diagrams, we find an incorrect relative factor of —2 between 
the infrared-singular and the finite terms in comparison to the gravity amplitude. 
This is presumably related to the lack of gauge invariance of the Wilson loop at the 
cusps. Motivated by these results, we then turn in Section 4 to consider a gauge where 
the cusp diagrams vanish, which we call the conformal gauge. We show that in this 
gauge the Wilson loop diagrams, where the propagator connects two non-adjacent 
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segments, precisely yield the full four-point Af = S supergravity amplitude, including 
finite and divergent terms, to all orders in the dimensional regularisation parameter 
e. This is in complete analogy to what happens in Af = 4 Yang-Mills in a similar 
gauge, as we show in Appendices lAl and [Bl 

Note added: After this work was completed, the preprint [75] appeared, which 
overlaps with Section [2] of this paper. 



2 MHV amplitudes in J\f = 8 supergravity and 
iterative structures 

In this section we start by briefly reviewing the expressions of the four-point MHV 
amplitude in A/" = 8 supergravity at one and two loops, and we then move on to study 
iterative structures at two loops. 



2.1 Background 



The form of the four-point MHV amplitude at L loops in maximal supergravity is 
very simple. It is given by the tree-level four-point MHV amplitude Ai^^^'^, times a 
scalar (helicity-blind) function, 

=MT"Mi^l (2.1) 

This amplitude was first calculated at one loop in [46] from the a' ^ limit of a 
string theory calculation, and later rederived in [47] using string-inspired techniques 
[48], as well as unitarity [49,50]. The infinite sequence of one-loop MHV amplitudes 
was obtained in [51]. Recently, the four- and five-point MHV amplitudes were also 
rederived in [52] using MHV diagrams. The two- and three-loop expressions were 
derived in [26], [35], respectively. 

At one loop, the function AA^^^ is simply given by a sum of three zero-mass box 
functions, 



where 



(1), , r dPi 1 

^4 [s,t) := y(^p(;_p^)2(^_p^_p^)2(;+p^)2 (2-3) 
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is a zero- mass box function with external, cyclically ordered null momenta pi, p2, ps 
and p4, which sum to zero. We set s := {pi +^2)^, t := {p2 +^3)^, u := {pi + PsY = 
—s — t, and D = 4 — 2e. Explicitly 



st 



cr 



TC 



(2.4) 



where cr := {47Ty-^r{l + e)r^{l - e)/r(l - 2e). Using (El, we can rewrite (E^D as 



M 



(1) 



Cr 



-u 



U 



1 — u log'' - — s log^ — — t log 



(2.5) 



The simplicity of (12.11) . where the tree- level amplitude factors out leaving a helicity- 
blind function of the particle momenta is clearly reminiscent of the structure for the 
infinite sequence of MHV scattering amplitudes in maximally super symmetric Yang- 
Mills. This motivates the search for 



a. an iterative structure in the higher-loop amplitude similar to that discovered 
in [1,2] for the A/" = 4 amplitude, and 



b. a derivation of the functions hA^^"^ using Wilson loops. 



The investigation of possible iterative structures of MHV amplitudes in A/" = 4 Yang- 
Mills was motivated by the known structure of the infrared divergences. This led 
BDS to propose in [2] the following conjecture for the all-loop MHV amplitude in 
A/" = 4 super Yang- Mills: 

00 00 
-Mn,YM := l+E«'-^SM(e) = exp[5^a^(/W(e)A^i^,^M(Le)+CW+E(^)(6))" , 

L=l L=l 

(2.6) 

where a = [g^N/{8TT^)]{ATre'^y. In ([2SD, f'^^Ke) = fj^^^ + f[^\ + /^^^e^ is a set of 
functions, one at each loop order, which make their appearance in the exponentiated 
all-loop expression for the infrared divergences in generic amplitudes in dimensional 
regularisation [8]. Specifically, /g'^^ = 7^''/4, where 7^- is the cusp anomalous di- 
mension (related to the anomalous dimension of twist-two operators of large spin). 
Importantly, the constants C*^^^ do not depend on the kinematics or on the number 
of particles n. The non-iterating contributions E'i^'* vanish as e — > and depend 
explicitly on n and the kinematics. 

As we have mentioned in the Introduction, the BDS proposal has been confirmed 
at two [1] and three loops [2] in the four-point case, at two loops for the five-point 
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amplitude [11], but the recent work of [12] shows a breakdown at two loops in the six- 
point case [12]. The infrared-singular part of A^„,ym is of course correctly reproduced 
by the infrared-divergent part of the right hand side of (12. 6p . 

In order to check (12. 6p . one takes the log and expands both sides in perturbation 
theory; for example, at two loops, one gets 



n,YM\ 



(1) 

n,YM\ 



(e))" = fHe)M^n!YMm + C(^^ + E(^K (2.7) 



We wish to follow the same path here for A/" = 8 supergravity, starting from the 
observation that in gravity the one-loop infrared divergences exponentiate [37]. In 
the four-point case, the leading infrared divergences are expected to resum to 



exp 



/K\ 2 2 / 

cr i^-j - i^s log(-s) + t log(-t) + u \og{-u] 



{21 



Notice the appearance of the invariant u = (pi + Ps)^, due to the lack of colour 
ordering. Moreover, in [51] it was shown that the tree-level soft and collinear splitting 
amplitudes in gravity are exact to all orders in perturbation theory. This is due to 
the fact that the coupling constant k is dimensionful, and it is always accompanied 
by a power of a kinematic invariant which vanishes in the limit considered [37,51]. 



We write the four-point MHV amplitude in TV" = 8 supergravity (stripped of the 
tree-level prefactor) aqj 



1 + 



4 



exp 



E 

.L=l 



where 



m 



(1) 



m 



(2) 



M 
M 



(1) 



(2) 



(1)^2 



(2.9) 



(2.10) 
(2.11) 



and so on. Motivated by (12. 6p and, specifically at two loops, by (12. 7p . we will calculate 
in the following section the difference appearing on the right hand side of (12. lip . 



Let us make a final comment before moving on to explore in detail iterative struc- 
tures at two loops. We observe that, unlike in the TV = 4 Yang-Mills case, the simplic- 
ity of (12. ip does not extend immediately beyond the four-particle the explicit 
results for the n-point amplitude of [51] show. It was shown in [51], using Af = 8 Ward 
identities, that the ratio , 2'^ , . . . , i~ , . . . , j~ , . . . , n~^) / {i j)^ is independent of 

■^Notice that in (|2.9p we absorb the appropriate power of k in the definition of A^4^'' and ■m'"/'\ 
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the positions z, j of the negative-hehcity gravitons, i.e. it is hehcity bhnd. This is 
similar to the Yang-Mills case [53], where J\f = 4 supersymmetric Ward identities 
allow one to move the position of the negative-helicity particle, and show that the 
corresponding ratio in Af = 4 Yang-Mills A^^(l+, 2+, . . . ,i~ , . . . . . . , (i j)^ 

is independent of i and j. In gravity however, this helicity-blind function is in general 
expressed as a sum of terms containing different spinor bracket valued coefficients. 
Two immediate consequences of this we would like to stress are that, firstly, it is 
not immediately clear what sort of iterative structures could be realised beyond four 
points; and, secondly, it is not obvious how a Wilson loop calculation could reproduce 
such terms (this situation somewhat parallels the problems one would encounter in 
attempting a derivation of non-MHV amplitudes in A/" = 4 super Yang-Mills from 
Wilson loops). For these reasons, in this paper we only concentrate on the four-point 
MHV scattering amplitudes. 



2.2 Iterative structure of the J\f = 8 MHV amplitude at two 



The previous discussion shows that, in searching for prospective iterative structures 
in the JV = 8 MHV amplitudes at two loops, it is meaningful to analyse the quantity 
(12. lip in supergravity, corresponding to the two-loop term in the expansion of the 
logarithm of the amplitude. We will carry out this computation in detail for the four- 
point MHV gravity amplitude described in the previous subsection. We observe that 
unlike the Yang-Mills ABDK conjecture [1] , but in agreement with Weinberg's result 
for gravity amplitudes [37], the one-loop infrared divergent terms of the amplitude 
exponentiate. More precisely, we will show that 



and calculate the function on the right hand side of (12.121) . 

The one-loop amplitude Ai^l^ is given in (12.21) . The two-loop amplitude was 
computed in [26], and is 



loops 




= finite. 



(2.12) 




( 



) 



4 



stu ■^( 



.,t) + .2jf'^(3,w) + .2jf)'NP(,,t) + ,2j(2).NP^^^^^^^y^j.^ 

(2.13) 
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Here I^^'^{s,t) and I^'^'''^'^ {s,t) are the planar and non-planar double box functions: 



K2),NP 



(27r)^ (27r)^ P {I - pif {I - Pi - P2f {I + kfk"^ {k - p^f {k - p^ - p^f ' 
d^l d^k 1 



{27t)d (27r)^ /2 (/ - p^y {I + ky{l + k + piY P {k - p3y {k-p3- piY 

(2.14) 



and in (12.131) we have to sum over the three cyclic permutations of the momenta p2, 
P3 and Pi (i.e. over the three cyclic permutations of s, t and u). 

The two-loop planar box function was first evaluated by Smirnov [54] (see also [2]) 
and the non-planar double-box function was evaluated by Tausk [55]. These expres- 
sions need to be evaluated in different analytic regions, due to the permutation of 
kinematic invariants: we fix s, t < but we will then need functions in which s or 
t are replaced by u = — s — t > 0, requiring a rather delicate procedure for analytic 
continuation. This procedure is outlined in Appendix ICl 

Smirnov's result for the planar double box integral (we use the form given in [2]) 
is given in terms of functions ^''^•''^(s, t) as 



Tf'^(.,t) = a,^ 
where := i (47r)'^~2r(l + e) and 



F(2).P(s,t) 



sH 



(2.15) 



-2e7 



r2(i + 6) 



Cj{-t/s) 



(2.16) 



with the coefficients Cj in (B.5) of [2]. This expression is valid in the region s,t < and 
we must carefully analytically continue into other regions as described in Appendix [Cl 

Tausk's expression [55] for the non-planar double box is given in terms of functions 
F(2)'NP(s,t) as 



.(2),NP 



(s, t) = a\ 



(2.17) 



The function ^^^^''^^(s, t) is given in [55] in all analytic regions (there it is called F^). 
Using the above results for the integrals, we arrive at the following expression for 
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the two-loop amplitude, 



M 



(2) 
4 



2 \ 2 

K a. 



SMF(2).P(s,t) + 2sMF(2).NP(s,t) 

+ t) + 2suF(2)'NP(w, t) + cyclic 



(2.18) 

Notice that the functions F^'^^'^{s,t) and F^'^^'^^ (s,t) always appear together in the 
combination F^"^^'^ + 2F^'^^'^^, although F^'^^'^{s,t) corresponds to the planar double 
box function (12.151) . whereas F^'^^'^^ {s,t) corresponds to one of the two terms in the 
non-planar double box function (12.171) . 



The one-loop amplitude (12.21) is expressed as a sum of zero-mass box functions 

F«(s,t)" 



lf\ where 



St 



and 



F«(.,t) 



-£7 



r 1 + 6 



-s 



i=-2 



Cj{-t/s) 



(2.19) 



(2.20) 



The coefficients cj are given in (B2) of [2]. Again this is valid for s,t < and we 
analytically continue to other regions. Together with (12.21) . this gives the following 
expression for the one-loop amplitude, 



u 



F^'\s,t)+tF^'\s,u) + sF^'\u,t) 



(2.21) 



On putting in the functions for all permutations - correctly defined in their re- 
spective analytic regions - into the formula for the amplitude (12.181) . we find that 
■Ma'^ — |(-^4^'*)^ is finite. This finite remainder is explicitly given in ( IC.Gp . As de- 
scribed in detail in Appendix [Cl this function can be considerably simplified to the 
following expression]^ 



M 



(2) 



K \4 

8^. 



[k{y) + k(l/y)] + [k(l -y) + k{l/{l - y))] 



+ e [k{y/{y -l)) + k{l- l/y)] + 0(e) , (2.22) 



''Notice that (|2.22p is somewhat formal, as there is no common region where aU the functions 
appearing are away from their branch cuts. The precise analytic continuations for the case s, t < 
are explained in detail in Appendix C, and the explicit, somewhat lengthier expression for the right 
hand side of (|2.22p valid in that region, is given in (jC.6p . 
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where 



+ i 



y + ^ - 45i,2(i/)L + - log^(l -y)+A S^M - 

Jtt \og\l -y)~ ^7r=^ log(l - - ALt: U2{y) + 47rLi3(i/) - 47rC(3) 



(2.23) 



where y = —s/t and L := log(s/t). Generahsed polylogarithms, including the Nielsen 
polylogarithms Sm,n which appear above, are discussed in [56]. 

After submitting this paper, we have compared our results to those of [75] , which 
contains a different form for the finite remainder fl2.22p . The two expressions are in 
fact in complete agreement. Specifically, one can rewrite fl2.22p as 



where 



st h{ — j + sth(^ ) + permutations 



^7rlog^(M;) - 87r5i,2(w) + 8 7rC(3) 



0{e) , 
(2.24) 

, (2.25) 



which after taking into account the different analytic regions considered (here we 
consider s,t < whereas the authors of [75] consider s,u < 0) is in precise agreement 
with the result of [75]. 



An interesting observation is that the functions appearing in the expression for the 
amplitude have uniform transcendentality. This is somewhat surprising - although 
the box function and the planar double box function have uniform transcendentality, 
the non-planar double box does not. Nevertheless, the combination of functions 
F^'^^'^^{s,t) + F^'^^'^^(u,t), which appears after summing over all permutation, does 
have uniform transcendentality. We notice that amplitudes in A/" = 1,4 supergravity 
do not have this property. This is explicitly shown by the calculations in [47] of 
the one-loop four-graviton MHV amplitudes, see Eq. (4.6) of that paper. Perhaps 
unexpectedly, the M = Q MHV amplitude is also maximally transcendental at one 
loop. It would be interesting to know if this property persists at higher loops in the 
perturbative expansion of the amplitudes in these theories. 



3 The one-loop Wilson loop calculation 

In this section we describe the one-loop calculation of the four-point MHV amplitude 
of gravitons from a Wilson loop. 
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The expression we are going to use is motivated by its apphcation in the eikonal 
approximation [43-45] to gravity [41,42], and it reads 



W[C] := (Pexp 



in (pdr h^i,{x{T))x^{T)x'' [t) 



c 



(3.1) 



where x^{t) parametrises the loop cjf| Note that the exponent in (13. ip can be rewrit- 
ten a^ 



'rf^xr'^^(x)V(a^) , (3-2) 
where, in the linearised approximation, the energy-momentum tensor is 

r'^"(x) := [dr x^'{t)x''{t)5^^\x-x{t)) . (3.3) 



The specific form of the contour C we choose is dictated by the graviton momenta 
Pi, ■ ■ ■ ,p4. In gravity there is no colour ordering - the amplitude (12.21) is a sum over 
the permutations (1234), (1243), (1324) of the four external gravitons. In order to 
match this from the Wilson loop side, we will therefore include the contribution of 
three Wilson loops with contours C1234, C1243, C1324, where Cijki is a contour made by 
joining the four graviton momenta pi, pj, pk, pi in this order. More precisely, the 
quantity we calculate at one loop will be 

W := W[Cu34] W[Cu43] W^[Cl324] • (3.4) 

Writing W[Cijki] := 1 + '}2T=i^^^K^ijki] = ^W^T=i'^ij'L the one-loop term of 
1^ is 

ly(l) = iy(l)[Ci234] + W^^'^[Cl243] + Vr(l)[Ci324] • (3.5) 

Before presenting the one-loop calculation, we would like to make a few prelimi- 
nary comments. 

1. One can check that the expression in (13. ip is not invariant under the gauge 
transformations 

h^u h^i, + df,^u + du^^ , (3.6) 

where ^'^(x) is an arbitrary vector field. Furthermore, it is easy to see that for contours 
composed of straight line segments joined at cusps such as those considered in this 
paper, the failure of gauge invariance is restricted to terms localised at the cusps. We 
think it is therefore not completely surprising that the infrared divergent parts of the 
Wilson loop will come out with an incorrect numerical prefactor from our calculation, 
compared to the finite parts, as we shall see below. 
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The same expression for the gravity Wilson loop has recently been used in [57]. 

uv- 



^In this section we set Z) = 4 — 2e 
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2. The expression (13.11) is not explicitly reparametrisation invariant, but it can 
be seen to arise from a reparametrisation invariant expression involving an einbein e, 
by writing the action of a free, massless particle as 

S ^ [ -T-\ i^^^^dtiv ■ 
J e(r) 

The energy momentum tensor resulting from this action is the one we use in our 
definition of the Wilson line in (13. ip . after gauge fixing e = 1. The equation of 
motion for the einbein just imposes the condition that the path of the particle is null. 
The contour of the Wilson loop we use is piecewise null so that no problems can arise 
from reparameterisation invariance away from the cusps. 

3. We note that the three contours appearing in (13.41) are obtained by permuting 
the external momenta, not the vertices. Due to the inherently non-planar character 
of gravity, one cannot consistently associate T-dual momenta to the external graviton 
momenta. For this reason, it is therefore unlikely that a version of dual conformal 
invariance might constrain the form of the amplitude here. 

4. A different expression for a gravity Wilson loop has been considered by 
Modanese [58,59], where the right hand side of (13.11) is replaced by 

(TrW(C)) , (3.7) 

where 



WJC) := Pexp 



c 



and is the Christoffel connection. The quantity TrW(C) has the advantage of 
being manifestly invariant under coordinate transformations [59]. The calculation 
of the one-loop correction to TrW(C) for a closed loop has been considered already 
in [59] , and the result is proportional to 

j^dx^^dy^^Tl^p{x)Tl^^{y)) . (3.9) 

We refer the reader to Appendix [D] for the details of the evaluation of (13. 9p in the 
linearised gravity approximation. The result is, dropping boundary terms, 

Idx^dy'' {Tlp{x)VtM) = c{D) Lx,dy^6'~''\x - y) , (3.10) 
Jc Jc 

where c{D) is a numerical constant which is finite as D ^ A. Parameterising the 
contour as x = x(cr), we can rewrite the right hand side of (I3.10p as 

c{D) (dr (da x^(r)x^(a) 5^^\x{t) - x{a)) . (3.11) 
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Some readers may notice that the divergent expression in (13.111) aheady appears in the 
lowest-order expansion of the Makeenko-Migdal loop equation [60,61] in Yang-Mills. 
An evaluation of (13. lip has been carried out in A/" = 4 super Yang-Mills in [62] for a 
cusped contour, by using a regularisation of the Dirac delta function which employs 
a cutoff of width a. Interestingly, the right hand side of (13.111) is then found to be 
proportional to l/a? times the one-loop cusp anomalous dimension^ We observe 
that, because of the delta function appearing in it, the expression in (13.111) receives 
contribution only from cusps and self-intersections present in the contour. The main 
point we would like to make here is that (13.111) does not reproduce (parts or all of) 
the A/" = 8 supergravity amplitude ( 12. 2p . Therefore, in the following we will work 
with the Wilson loop defined for a polygonal contour as in (13. ip . 

We now proceed to describe the calculation. We work in the de Bonder gauge, 
where the propagator is given by 

1 / 2 \ 

(^^Jl/X2(^)^I/l 1^2(0)) ~ lj^\JltJ.lVlVtl2f2 ~^ VfJ.lV2VtJ.2Vl ~ JJ _ 2^Ml/^2^l'l!^2 1 (3.12) 

where 

■= -^r(--l) (3.13) 



w ) 



The gravity calculation is very similar to the one-loop calculation performed in [17,18] 
for the one-loop Wilson loop in maximally supersymmetric Yang-Mills theory. As in 
that case, three different classes of diagrams contribute at one loopj^ In the first one, 
a graviton stretches between points belonging to the same segment. As in the Yang- 
Mills calculation, these diagrams give a vanishing contribution since the momenta of 
the gravitons are null. In the second class of diagrams, a graviton stretches between 
two adjacent segments meeting at a cusp. In the Yang-Mills case, such diagrams lead 
to ultraviolet divergences [64-71]. As in the Yang-Mills Wilson loop case [17], these 
divergences are associated with infrared divergences of the amplitude by identifying 

£ttv — £■ 



We will now see how in our gravity calculation, these divergences are still present 
but will be softened (from 1/euv^ to 1/euv) after taking into account the sum over 
the contributions of the three Wilson loops. 

A typical diagram in the second class is pictured in Figure [H There one has 

^In [63] this result was extended to two loops, and conjectured to hold to any loop order in 
perturbation theory. 

^For a Wilson loop bounded by gravitons, only gravitons can be exchanged to one-loop order. 
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Figure 1: A one-loop correction to the Wilson loop bounded by momenta pi, ■ ■ ■ ,p4, 
where a graviton is exchanged between two lightlike momenta meeting at a cusp. Dia- 
grams in this class generate infrared- divergent contributions to the four-point ampli- 
tude which, after summing over the appropriate permutations give rise to f l3.16p . 



xi{ti) — X2{t2) = pi{l — Ti) + P2T2- The cusp diagram gives 



■(m/i 



-(zK/i'uv)2 



r(i - 6„v) 

47r2-euv 

r(i - 6,,v) 

47r2~<^uv 



dTidT2 — ^ 

[-{piTi+P2r2) Y- 

1 r-s)i+^uv- 



(3.14) 



Notice that we need to choose euv > in order to regulate the divergence in fl3.14p . 
Furthermore the scale used in the Wilson loop calculation is related to the scale used 
to regulate the amplitudes fi a.s ft = (c/i)~^ (the precise coefficient c in front of fi can 
be reabsorbed into an appropriate redefinition of the coupling constant). 



Summing this over the four cusps of the first Wilson loop, one get^ 

c(euv) 



^tuv 



-t) 



1+euv 



Adding the contributions of the two other Wilson loops, we get 



(3.15) 



(3.16) 



Upon expanding this expression in euv, the cancellation of the 1/euv^ pole becomes 
manifest (after using s -\- 1 + u = 0), and (13.161) becomes, up to terms vanishing as 



^We set c(euv) = (K/i'uv)2r(i - euv)/(47r2-^uv). 
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euv 0, 



c(euv) — (slog(-s)+tlog(-t)+Mlog( 



-u 



s log^(-s) +t log^(-t) +u log^ 



-u 



(3.17) 

We recognise that this expression is the infrared- divergent part of the four-point MHV 
gravity amphtude (12.21) . We notice however that, after summing over the appropriate 
permutations as in (13. 4p . one finds that these infrared- divergent terms have an extra 
factor of —2 compared to the finite parts, to be calculated below. We believe this 
mismatch is not unexpected, given that the failure of gauge invariance of (13. ip occurs 
at the cusps0 

X3 



X2 




X4 



Figure 2: Diagrams in this class, where a graviton stretches between two non-adjacent 
edges of the loop, are finite, and give in the four-point case a contribution equal to 
the finite part of the zero-mass box function F^^\s^ t) multiplied by u. 

We now move on to the last class of diagrams, where a graviton is exchanged 
between two non-adjacent edges with momenta p and g; one such example is depicted 
in Figure [2l In the Yang-Mills case these diagrams were found to be in one-to- 
one correspondence with the finite part of the two-mass easy box functions with 
massless legs p and q. We will show now that (13. ip leads exactly to the same kind of 
correspondence with the finite part of the one-loop four-graviton amplitude. 

Indeed, the one-loop diagram in Figure [2] is equal to 

c(e,v) fdndr^ — ^^^^ . (3.18) 

Jo [-(Pl(l - Ti) +P2 +P37-2) ]^-^uv 



A factor of 2 could be explained because we are effectively double-counting the cusps in summing 
over the permutations, however at the moment we are unable to explain the relative minus sign. 
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This integral is finite in four dimensions, and gives 



(3.19) 



Summing over the two possible pairs of non- adjacent segments and including the 
contributions of the two other Wilson loop configurations, we get exactly the finite 
part of the one-loop MHV amplitude in A/" = 8 supergravity (12.51) up to the tree-level 
amplitude0 



4 Calculation in the conformal gauge 

The gravity Wilson loop defined above, unlike the Yang-Mills Wilson loop, is gauge 
dependent. It turns out that one can define a gauge in both cases in which the cusp 
diagrams vanish completely. We call these "conformal" gauges!^ In the Yang-Mills 
Wilson loop one obtains the same answer in either gauge, but in the gravity Wilson 
loop the conformal gauge appears to be the unique gauge which gives the amplitude, 
both infrared- divergent and finite pieces correctly, to all orders in e. 



4.1 Gravity propagator in general gauges 

We first need to define a general class of gauges in the gravity case. To do this, we 
consider the free Lagrangian of linearised gravity: 



which can be easily checked to be invariant with respect to the gauge transformation 
5h^u = 2d(^^^y) . We then add a gauge fixing term of the following form: 



^^A Wilson loop calculation clearly cannot produce any dependence on helicities and/or spinor 
brackets. Incidentally, we also observe that in Yang-Mills, a Wilson loop calculation cannot produce 
any parity-odd terms such as those appearing in the five- and six-point two-loop MHV amplitudes. 

^^This name is motivated by the fact that, in the Yang-Mills case, the D-dimensional propagator 
turns out to be proportional to the inversion tensor J^j/(x) :— rj^v — 2xp,Xi,lx^. The Yang-Mills 
conformal propagator is described in Appendix \^ where we show that it can be obtained from 
a Feynman-'t Hooft gauge- fixing term with a specific coefficient. In Appendix IB] we perform the 
calculation of the n-point polygonal Wilson loop. The outcome of this calculation is that cusp 
diagrams in the conformal gauge vanish, and the TV = 4 amplitude is obtained from summing over 
diagrams where a gluon connects non-adjacent edges. In this case, each such diagram is in one-to-one 
correspondence with a complete two-mass easy box function. 



(4.1) 
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which is de Donder-hke, but with an arbitrary free parameter a. We will call this the 
a-gauge. 



In momentum space, the corresponding gauge-fixed Lagrangian has the form 
(l/2)/i'^'^i^^,,^v'/i'^'"', where 



K^u,fM'u'{k) = k^r]^,(^^r]^y - 2ki^^r]^)^^>k^/) - k'^rj^yrj^i^i + r]^^k^>k^' + v^'fc^fc, 
— a 



k{f,ri^){u'k^>) - ^{i]^ukf,'k„> + r]^'u'kf,ku) + ^k^VfiuVfi'w ■ (4.3) 



Now we define the propagator D^^y^i to be the inverse of K^i,y^i, i.e. 

i^M-,^'.'/^'^'"' = S^^d'l^ . (4.4) 

By writing down the most general Lorentz covariant terms which have the correct 
index symmetries and have mass dimension equal to -2, we see that D^^^^iyi must 
take the form 

c d 
+ -j^iVf^uk^^'k^' + Vfi'u'k^k^) + -j^kf,kykf,iky> . (4.5) 

Then (14.41) gives a set of equations for the free parameters which have the unique 
solution (for D ^ 2), a = -(4 + 2a)/a, 6 = -!/(£> -2), c = ci = 0. Thus, the 
propagator corresponding to the a-gauge defined by the gauge-fixing term (14. 2 p is 
given by 



/,A W 1 \ 4 + 2a 1 , , ,^ ^, 

^''^'(^^ = -j^yl^^'i^^Vu)u' - jj^Vf^uVt^'w ) — -j^hf.Vu)iu'k^,>) . (4.6) 

Notice that (14.61) reproduces the standard de Bonder propagator for a = —2. 



4.2 Propagator in position space in the a-gauge 



We now perform the Fourier transform to position space. The Fourier transform of 
1 / k'^^ has the fornix 

= c{D, A)(-x2)^-^/2 , (4.7) 

where 

^ ' ' r(A + 1 -D/2)r(A) ^ ' 



^•^More details can be found in Appendix [XI 
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By differentiating twice with respect to x and setting A = 2 we find that the Fourier 
transform of k^ku/k'^ is 



2c(D,2)e, 



+ . . , (1 - euv) 



2')l— SUV ( x^)^^^'^^ 



(4.9) 



Using this we take the Fourier transform of fl4.6p . and obtain the propagator in 
position space: 



V(m^^)^' 4D, 1) 1 



-X- 



1 



-X 



2\2-evv 



(4.10) 



where 



A = c(D,l) + 2aeuvc(D,2) = 4a euv(l - Cw) c(Z^, 2) , (4.11) 



and 



4 + 2a 



a 



(4.12) 



4.3 The conformal gauge 

By direct analogy with the Yang-Mills case, discussed in Appendix [Bl where we 
show that in the "conformal" gauge the cusp diagrams vanish, we define the gravity 
conformal gauge to be the gauge in which the cusp diagrams vanish. We show in this 
section that this particular gauge can be obtained from an a-gauge fixing term as 
defined in the previous section for an appropriate value of the parameter a. 

To begin with, consider the cusp defined by momenta p, q and then let x = pa+qr. 
Then the term appearing in the cusp at one loop is 

P^p'^M-X-Ka;)// = {-x'y^^'^pqfariB - 2A) . (4.13) 

Therefore, the cusp diagrams vanish for B = 2A. One can quickly check that this 
implies a = —c{D, l)/(2euv^c(D, 2)) = 4/(D — 4). The result is the propagator in the 
conformal gauge: 



D^iu,fi'u'{x) = c(L),l)^'^^ 



e 



uv 



1 



(4.14) 



which requires 

a = -2(L) - 4)/(L) - 2) . (4.15) 
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4.4 Gravity Wilson loop in the conformal gauge 



We now proceed to calculate the gravity Wilson loop in this conformal gauge. We 
have shown that the cusp diagrams are equal to zero in this gauge, therefore we 
need only calculate the "finite" diagrams (which are now no longer finite). Consider 
the Wilson loop with edges Pi,P2,P3,P4 (in that order) and the graviton stretching 
between sides 1 and 3. Then we have x = api + rpa + p2 and = sa + tr + uar. 
The contribution of this diagram is then 

/ dadTpiVD^.u,f.'u'{x)p3^''p/ (4.16) 
Jo 

= c{D,l)^^ — / dadr- 



euv 4 Jo {-{so + tT + uaT)y- 

1 u 



c{D, 1) - {-sy'^^2Fi{l, euv, 1 + euv, 1 + 7) - (-t)^"^2Fi(l, Cuv, 1 + Cuv, 1 + -) 

6uv 4 L t S - 

We see that we obtain the complete (infrared-divergent as well as finite pieces) two- 
mass easy box function to all orders in euv Adding the other diagram (which gives 
the same result) and then summing over the remaining permutations as described 
above, gives the correct one loop A/" = 8 supergravity amplitude fl2.2ip . 

Despite this encouraging result, we should remember that our starting expression 
for the Wilson loop (13.11) was not gauge invariant. It would be important to remedy 
this gauge non-invariance, which is localised at the positions of the cusps, by an 
appropriate subtraction procedure. Furthermore, it would be interesting to study 
infrared divergences, as well as the derivation of finite parts of gravity amplitudes at 
higher loops using the Wilson loop proposed in (13.41) . 
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A The conformal propagator in Yang-Mills 



In this section we briefly outline the construction of the conformal propagator. It is 
defined to be proportional to the inversion tensor 



By using 



1 



TT 



2 /D 



X^ 



(27r; 



D 



ipx P/J-Pi^ 
p4 



4-^2 

vr 2 ^/D 



-x^ + le) 2 



— -1 



8 V 2 



r]f,y - {D - 2)x^Xy/x'^ 



-x^ 



le) 2 



it is easy to see that the following combination has the desired property: 



d^p 



ipx V/J-i^ _|_ 

(27r)^ p2 



where we define the conformal propagator 



d^p 
(2^' 



^ipx 



PpPu 
p4 



r 



D 



D-A 4 



-x^ + ie) 2 ■ 



X ^Xi, 



x^ 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



Thus, the expression f lA.4|) is obtained by choosing a Feynman-'t Hooft gauge-fixing 
term (a/2) J d^a; (9^A^)^ for the particular choice of a = {D — A)/ D. The vanishing 
of this gauge-fixing term in D = 4 dimensions is reflected in the presence of a factor of 
1/(D — 4) in (lA.4p . which makes this propagator not well defined in four dimensions. 



B The Yang-Mills Wilson loop with the conformal 
propagator 

As a simple but illuminating application of the above conformal propagator, we would 
like to outline the calculation of the Yang-Mills Wilson loop with a contour made of 
n lightlike segments performed in [18]. Of course, the usual expression of the Wilson 
loop in Yang-Mills is gauge invariant, hence evaluating it in any gauge leads to the 
same result. The use of this gauge leads however to a recombination of terms, where 
the cusp diagrams vanish!^ Consider for instance the cusped contour depicted in 
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Figure 3: A one-loop correction for a cusped contour. We show in the text that, when 
evaluated in the conformal gauge, the result of this diagram vanishes. 

Figure [31 Using the conformal propagator, and Xp^(ri) — Xp2{T2) = — ri) + P2T2, 
we see that the one-loop correction to the cusp is given by an expression proportional 
to 

IIU _ 2 [Pl(l-n)+P2T2]^[pi(l-ri)+p2T2]'^ 
J J 2(pip2)(l-ri)T2 ^ N 

dTidT2 pii,p2u r^TT rjz ^ 1^/2-1 ' 

[-2{pip2){l - TijrsJ^/^ ^ 

which vanishes. 

We now move on to consider diagrams where a gluon is exchanged between non- 
adjacent segments, such as that in Figure |H In [18] it was shown that this diagram 
is equal to the finite part of a two-mass easy box function. In the conformal gauge, 
a simple calculation shows that it is equal tc0 

1/ „o dr-i dro , 

where 

D{n,T2) := {x,{n) - {x,{r2)f (B.3) 

= P^ + is- p2)(l - n) - (t - P^)t2 - u{l - n)T2 , 

where we used 2{pP) = s ~ P^, 2(gP) = t- P'^, and s + t + u = P^ + Q^. We have 
also introduced 



""^^The usual infrared-divergent terms are produced by diagrams which, in the Feynman gauge 
calculation of [18], were finite. 
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In the following we set e = — 
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Figure 4: A one-loop diagram where a gluon connects two non-adjacent segments. In 
the Feynman gauge employed in [18], the result of this diagram is equal to the finite 
part of a two-mass easy box function F^^^{p, q, P, Q), where p and q are the massless 
legs of the two-mass easy box, and correspond to the segments which are connected by 
the gluon. In the conformal gauge, this diagram is equal to the full box function. The 
diagram depends on the other gluon momenta only through the combinations P and 
Q. In this example, P = pz + Pi, Q = Pe + Pr + Pi- 



In [18] it was found that 



dri dT2 



[-D(ri,r2)]2+^ P2 + g2_,_i 



where 





1 




" ~^ 




a \ 


(r 


-aPV 




a y 


(r 


— as) 



(B.5) 



(B.6) 



e, e, 1 + e 



1 



' 1 - aP2 
1 

1 — as 



1 



a 



l-ag2 

2F1 ( e, e, 1 + e, 



1 — at 



1 — at 



where we have introduced 



a :- 



P^ + Q^-s-t 



(B.7) 



P2g2 _ 

Notice that in fIB.Sp the function JF appears with argument e + 1. After a moderate 
use of hypergeometric identities, we find that the one-loop correction in flB.2p is equal 
to 

1 rn -^ . . . 

(B.8) 



ir(l+e) ^2me(^^^^p2^g2^ 



2 47r2+^ 
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where F'^'^^{s,t, P'^,Q'^) is the all-orders in e expression of the two- mass easy box 
function derived in [72] 



.p2 



/i 



2 



2Fi(l,-e,l-e,aP2) 



2F1 (1, -e, 1 - e,as) + 



2Fi(l,-e,l-e,at) 



2F1 (l,-e,l-e,aQ2) 



(B.9) 



Summing over all possible gluon contractions in the Wilson loop, one finds complete 
agreement with the result derived in [18] for the same Wilson loop, as anticipated. 



C Analytic continuation of two-loop box functions 



In Section[2]the one and two loop amplitudes are given in terms of functions t), 
F^'^^'^^(s,t) and F^^\s,t). In Yang-Mills, colour ordering means that we need to de- 
fine the functions explicitly in only one analytic regime. In gravity however, we must 
sum over permutations of the kinematic invariants. Even if we fix the kinematic 
regime to be s, t < we must also consider for example F{s, u), and the second argu- 
ment of this function will be greater than zero (recall that u = —s — t). There will be 
three different kinematic regimes of interest and, following Tausk [55], we label them 
in the following way: 

{Fi{s,t) t,u<Q 
F2(s,t) s,u<0 (C.l) 
F3(s,t) s,t<0 . 



Tausk gives explicit formulae for the non-planar box function in all three regions, 
but it is nevertheless useful to know how to obtain the function in any region from 
its manifestation in a particular region. The Mathematica package HPL [74] is very 
useful for this. 

We will sketch the procedure below. Let us begin by considering the analytic 
continuation from region 1 to 2. In general, functions in this region take the following 
form: 

F,{s,t) = f{\og{s),\og{-t),\og{-u),Hs,i{-t/s)) . (C.2) 

Here H^^i{z) is a harmonic polylogarithm where a represents a string of zeros or ones. 
Note that at two loops we need not use harmonic polylogarithms as they can all be 

i^Omitting a factor of cr = r(l + e)r'^{l - e)/(47r)2-^ compared to [72]. 
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re-expressed in terms of Nielsen polylogarithms. On the other hand, at higher loops 
harmonic polylogarithms will appear which cannot be so expressed; it is nevertheless 
useful to use harmonic polylogarithms even here (see [56, 74] for more details on 
harmonic polylogarithms) . Such a harmonic polylogarithm is analytic everywhere on 
the complex plane except for a branch cut on the real axis for z > 1. Note that the 
arguments of all the (poly)logarithm functions in (lC.2p lie away from the branch cut. 

Now the function continued to region 2 takes the following form: 

F2(s,t) = /(log(-s) +z7r,log(t) -m,log(-M),i/5,i(-t/s)) . (C.3) 

We have analytically continued the logs appropriately, however the argument of the 
HPL functions now lies on the branch cut in region 2 {—t/s = 1+u/s > 1). We use the 
HPL package to transform away from the branch cut. Specifically putting —t/s = 1/y 
the command 'HPL Convert To Simpler Argument' will rewrite this in terms of HPLs 
with the argument y = —s/t which lies off the branch cut (one must also use the 
command 'HPLReduceToMinimalSet' to write the functions in a standard form). 

If we wish to obtain the formula in region 3 from that in region 1 we immediately 
have a problem. The argument of our HPL functions is —t/s which is not on a branch 
cut for either region. However, close examination shows that as we pass smoothly 
from region 1 to region 3, we must first pass along the branch cut - for example we 
must pass through the point s = 0, i.e. —t/s = oo. The HPL programme will not 
take this into account and the naive analytic continuation gives the wrong result. So 
it is better to first perform a transformation y ^ 1 — y on the HPLs in Fi{s, t) to find 
a new expression for -Fi(s, t) in terms of HPLs with argument 1 +t/ s = —u/s, i.e. 

F,{s,t) = g{\og{s),log{-t),\og{-u),H,,,{-u/s)) . (C.4) 

Then in region 3 we find —u/s > 0, and hence we are on the branch cut and we can 
proceed as before. We analytically continue as follows, 

Fs{s, t) = g(\og{-s) + ZTi, log(-t), log(-M) - in, Hu,i{-u/s)) . (C.5) 

Now use the HPL programme to transform back off the cut using the transformation 
y ^ 1/y yielding HPLs with argument —s/u. 

Now we have found the functions in all three analytic regions, and we can trans- 
form the arguments to obtain all the different permutations entering in the two-loop 
amplitude fl2.13p . For example F{s,t) = F3{s,t) since we are in the region s,t < 0, 
but F{u,t) = Fi{u,t) since the first argument is positive etc. 

At this point, after summing all contributions, the two-loop amplitude will be a 
linear combination of harmonic polylogarithms with different arguments. We there- 
fore use the HPL programme again to transform them all to the same argument. 
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ensuring that we never land on a branch cut in so doing. For example, for harmonic 
polylogarithms of the form Hg^i{x) (i.e. where the defining string of numbers ends in 
a '1') we restrict ourselves to transformations of the form y ^ 1 — y and y ^ y/{y — l) 
which the HPL program performs assuming we are away from the branch cut. 



Using the above techniques we obtain the following form for the two-loop finite 



remainder — h{M^^^)'^: 



2 



(1)a2 



2 N 2 

K a. 



's'f^'\y)+ef^'\y)+u'f^-\y) 



where 

f^'\y) =y - ^ log(l - y)L' - \og\l - y)L' + n'L' + ^ V(l - y)L 

- log(l - y)L + SS^^L - 47r2Li2(y) + 8S,M - SS^M 



+ i 



27rL3 



+ 27rlog(l - y)L^ - 27rlog2(l - y)L + 8TTU2{y)L 



+ 



47r^L 



-87rLi3(y) + 87r5i,2(y)] , 



(C.6) 
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(C.7) 



I log(l - y)L' + log^(l - y)L^ + 4U2{y)L^ - n^L^ - \ \oi{\ - y)L 



+ 47r' log(l -y)L- %lAj,{y)L + ^-K^lA^iy) + 8Li4(y) - 8>Si,3(y) + 



TT 



+ i 



2'kV' 



- 27rlog(l - y)L^ + 27rlog^(l - y)L 



-87r5i,2(?/) + 87rC(3) 



(C. 



f(-){y) ^1 log4(i -y)- ^Llog^(l -y)+L' \og\l -y)- ^L' log(l - y) - AL'U^iy) 



TV 



8LLk{y) - 8LU{y) - 8LSM + 8^2,2(1/) - y + L^tt^ 



+ i 



27rlog(l - y)L^ + 27rlog^(l - y)L - 87rLi2(y)L 



+ ^ - ^vr V(l -y)- ^nHog{l - y) + 8nUs{y) - 8nC{3) 
where y — —s/t and L :— \og{—y). 



(C.9) 



Since the amplitude is invariant under crossing symmetry (arbitrary permutations 
of the momenta or equivalently arbitrary permutations of s, t, u) we must have 

= /W(l/z/) = /(^)(1 -y)= f^\y/{y - 1)) , (C.IO) 
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which one can indeed verify as long as one takes suitable care over the analytic 
continuation in the manner outlined above. 



Simplifying slightly f^^\y) by writing it as k{y) + k{l/y) we obtain the form of 
the amplitude given in fl2.22p . 



D Derivation of ( l3.1Ulj 



In this Appendix we derive (13.101) from (13. 9p in the linearised gravity approximation. 
Upon expanding the metric about fiat space, gfj,u{x) = rj^u + K'h^u^x), one finds that 
(13. 9p is equal to 



c 



c 



To perform the calculation in fID.ip we choose the de Bonder gauge, where the propa- 



gator in D = 4 — 2euv dimensions is given by (I3.12p . Boundary terms can be dropped 
as the contour is a closed loop. Doing this, one easily finds tha10 

jdx^dy'' {Vlp{x)VlM) = c{D) jdx^dy^^U,^{x-y) 

= c{D)j^ dx^dy^5^''\x-y) , (D.2) 
where c{D) is a numerical constant, finite as ^ 4. This is the result quoted in 



^^In [59], terms such as those appearing on the right hand side of (|D.2p are referred to as "ultra- 
local". 
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